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Abstract 



Tenuous dust clouds of Jupiter's Galilean moons lo, Europa, Ganymede and Callisto have been de- 
tected with the in-situ dust detector on board the Galileo spacecraft. The majority of the dust particles 
have been sensed at altitudes below five radii of these lunar-sized satellites. Wc identify the particles 
in the dust clouds surrounding the moons by their impact direction, impact velocity, and mass distri- 
bution. Average particle sizes arc between 0.5 and 1 /y,m, just above the detector threshold, indicating 
a size distribution with decreasing numbers towards bigger particles. Our results imply that the par- 
ticles have been kicked up by hypervelocity impacts of micrometeoroids onto the satellites' surfaces. 
The measured radial dust density profiles are consistent with predictions by dynamical modeling for 
satellite ejecta produced by interplanetary impactors (Krivov et al., Planet. Sp. Sci., 2003, 51, 251- 
269), assuming yield, mass and velocity distributions of the ejecta from laboratory measurements. A 
comparison of all four Galilean moons (data for Ganymede published earlier; Kriiger et al., Planet. 
Sp. Sci., 2000, 48, 1457-1471) shows that the dust clouds of the three outer Galilean moons have 
very similar properties and are in good agreement with the model predictions for solid ice-silicate 
surfaces. The dust density in the vicinity of lo, however, is more than an order of magnitude lower 
than expected from theory. This may be due to a softer, fluffier surface of lo (volcanic deposits) as 
compared to the other moons. The log-log slope of the dust number density in the clouds vs. distance 
from the satellite center ranges between -1.6 and —2.8. Appreciable variations of number densities 
obtained from individual flybys with varying geometry, especially at Callisto, are found. These might 
be indicative of leading-trailing asymmetries of the clouds due to the motion of the moons with respect 
to the field of impactors. 
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1 Introduction 



All celestial bodies without atmospheres are permanently exposed to bombardment by hypervelocity 
micrometeoroids which knock-off secondary ejecta dust particles from the surfaces of these bodies. Im- 
pact ejection of dust particl es has been suggested as the main process for ma i ntaining dusty planetary 



Burns et al 



rings like the Jovi an rings ()Morfill et all 1198(1: iHoranvi and Cravena. Il99(i: lOckert-Bell et 



19991 ^. Saturn's E ring (IHoranvi et al 
well as putative dust belts of Mars Jsoteil. Il971 



19991) and 



luto (IThiessenhusen et al 



spacecraft (|Grun et all Il99^ a dust cloud form ed by impact ejecta particl es was for the first time 



1992 



_al 



Hamilton. D. P. and Burns. J. A 



Krivov and Hamilto 



nl. ll997l : 



199fl 



19941'! as 



Krivov and .Turewicz , 



20021 ). With the in-situ dust detector on board the Galileo 



detected surrounding Jupiter's moon Ganymede (|Kruger et al 



1999^ 



i 



20001) . The dust cloud was by 



far too tenuous to be detectable with remote sensing techniques. Particles belonging to the Ganymede 
dust cloud were identified by their impact direction and impact speed and their sizes were mostly 
below 1 fim. The spatial distribution of the grains as well as their size distribution were in agreement 
with model predictions based on the impact ejection mechanism. 

The Galileo dust measurements can be treated as a natural impact experiment leading to the detec- 
tion of the ejecta of hypervelocity impacts in space. They can give more insight into the process of 
hypervelocity dust ejection, for which the laboratory experiments on Earth still do not yield a com- 
prehensive picture. The measurements of the dust cloud at Ganymede stimulated the development 
of analytical models for impact-gener a.ted circumsa t ellite dust clouds not only for the Galilean moons 



but also for the Saturnian satellites (Krivov et al 



2003; 



Sremcevic et al, 



20031 ) . Thi s is especially 



impo rtant for the dust measurements to be collected at Saturn with the dust instrument (|Srama et al 
2OO2I ) onboard the Cassini spacecraft beginning in 2004. 



Since December 1995 Galileo has been on a bound orbit about Jupiter. The spacecraft had a total 
of 32 targeted close fiybys at all four Galilean moons: 7 encounters with lo, 11 with Europa, 6 
with Ganymede and 8 with Callisto. During many of the encounters between 1995 and early 1999, the 
impact rate o f dust grains showed a sharp peak within about hal f an hour centered at closest approach 

1999a|) . These peaks indicated the existence 



to the moon (Griin et al 



1996 



1997 



199 



Kriiger et al 



of dust concentrations not only at Ganymede but also in the close vicinities of lo, Europa and Callisto. 
During the fiybys at the Galilean moons after mid-1999, the spacecraft orientation prevented the 
detection of dust particles close to the satellites. In November 2002, Galileo had th e only opportunity 
to in -situ measure dust in the close vicinity of a fifth Jovian moon: Amalthea ( Kriiger and Griin . 
20oi). 



In addition to dust clouds surrounding the Galilean moons, at least three other populations of dust were 
detected by Galileo in the Jovian system (IGriin et allll998l ). Streams of 10-nanometer dust particles 
were detected throughout t he Jovian magnetosphere and were rec ognizable even in interplanetary space 



out to 2 AU from Jup it er (IGriin et al 



( Horanvi et al 



1993 



alibi: iGraps 



et al 



1993: 



200c 



Griin et al 



Kriiger et al 



1994 ). These dust grains originate from lo 



2003a|), their ultimate source probably being 
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the most powerful of lo's volcanic plumes. Bigger, micrometer-sized particles form a tenuous dust ring 
between the Galilean moons and further away from Jupiter. Many of these particles orbit Jupiter 



on prograde orbits whereas a population on re trograde orbits exists as well (|Colwell et al 



Thiessenhusen et al.l . 1200(1 : iKrivov et al 



. l2002afl 



1998b 



An overview of the Ga. 



Jupiter including the dust instrument itself can be found in 



Kriieer et al 



ileo dust measurements at 



(2003^) 



A detailed analysis of the dust grains detected at Ga lileo's four Ganymede flybys in 1996 and 1997 



has been published earlier (jKriiger et al 



2 



2OOOI ). showing that this Jovian moon is surrounded 



by a dust cloud formed by impact ejecta. Here, we analyze the dust impacts detected close to lo, 
Europa and Callisto, and compare our results with the measurements at Ganymede. Relevant physical 
properties of these moons are summarized in Table In Section [21 we give a brief overview of the 
most important aspects of the Galileo dust instrument, the Galileo spacecraft and the procedure to 
identify impacts of ejecta cloud particles in the Galileo dust data set. In Section |31 we analyze the 
dust detections at lo, Europa and Callisto. In Section^ we compare the properties of the dust clouds 
of all four Galilean moons and check them against modeling. Section El lists our conclusions. 



Insert Table HI 



2 Dust impact detection 



2.1 Galileo dust instrument 



The Galileo in-situ dust measur ements at Ganvmede p rovided the first in-situ detection of an impact- 
Processing of these measurements has been 



generated dust cloud in space (iKriiger et al 
described in detail by 



Kriiger et al 



i2OO0l ). 



1999 



Here we apply the same analysis techniques to the dust 
measurements obtained in the close vicinities of Europa, Callisto and lo. Descriptions of the dust 
instrument, Galileo spacecraft, data transmission etc. have been published in previous papers. In 
what follows we recall only the most important aspects and give references to earlier publications 
where necessary. 

Galileo is a dual-spinning spacecraft, with an antenna that points antiparallel to the positive spacecraft 
spin axis. During most of Galileo's orbital mission about Jupiter, the antenna pointed towards Earth. 
The Galileo Dust De tector System (DD S), like its twin on-board Ulysses, is a multi-coincidence impact 



ionization detector (Griin et al 



19921 ) which measures submicrometer- and micrometer-sized dust 



impacts onto the detector target. The dust instrument is mounted on the spinning section of Galileo 
and its sensor axis is offset by an angle of 60° from the positive spin axis (Fig. Thus, during 
one spin revolution of the spacecraft, the detector scans the entire anti-Earth hemisphere, whereas 
particles approaching from the Earth-ward direction remain undetectable. 

For each dust grain hitting the detector target, three independent measurements of the impact-created 
plasma cloud are used to derive the impact speed v and the mass m of the particle. The charge Qi 
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relea sed upon impact onto the target is described by the relation (|Goller and GriinLll989l : lGrun et al. . 

mi 



Qiocm-v^-^. (1) 

The calibrated speed range of the dust instrument is 2 to 70kms^^. The coincidence times of the 
three charge signals together with the charges themselves are used to classify each impact into one 
of four categories. Class 3 imp acts have three charge signals, two are required for cl ass 2 and class 1 
events, and only one for class oC 



iguhill99.4 



Griin et al 



1995; 



Kriiger et al 



1999H ). Class 3 signals, 



our highest class, are real dust impacts a nd class eve r its are noise. Class 1 and class 2 events are 
true dust impacts in interplanetary space (jBaguhl et al.l . 119931 : iKriiger et al.l . ll999bl ). However, in the 
Jovian system, within about 15Rj distance from Jupiter, energetic particles from the Jovian plasma 
environment cause an enhanced noise rate in class 2 and the lower quality classes. By analysing the 
properties of the lo stream parti cles and comparing t hem with the noise events, the noise could be 
eliminated from the class 2 data ( Kriiger et al.l . |l999g) . All class and class 1 events detected in the 
Jovian environment are usually classified as noise. 



Our noise identification scheme (jKriiger et al.l . ll999d ). however, was derived for the Jovian dust stream 
particles and, hence, its applicability to other populations of dust had to be verified. Since Europa 
orbits Jupiter within the region where the high noise rates occurred, a slightly modified scheme has 
been developed for the ejecta grains detected in the close vicinity of Europa (see 



Kriiger et al 



2001 



their Table 4). It will be applied in this paper to remove noise events from the data sets obtained at the 
flybys at Europa and lo. For the Callisto data, no noise removal is necessary because Callisto orbits 
Jupiter outside the region where the high nois e rates occurred. N oise removal was also not necessary 
in the earlier analysis of the Ganymede data ( Kriiger et al. . 200o[ ). It has to be noted that the noise 



removal technique uses statistical arguments and is applicable to large data sets only. Individual dust 
impacts may be erroneously classified as noise and vice versa. 

Galileo has a very low data transmission capability because its high-gain antenna did not open com- 
pletely. For the dust measurements this means that the full set of parameters measured during a 
dust particle impact (spacecraft rotation angle, impact charges, charge rise times, etc.) could only be 
transmitted to Earth for a limited number of impact events. During the close satellite flybys at the 
Galilean moons these limits were between one event per minute (record mode) and one event per 21 
minutes (real time science mode) (jKriiger et alJ ■ l2001^ . When event rates (i.e. dust impacts plus noise) 
exceeded these numbers, the full set of parameters was transmitted to Earth for only a f raction of all 



detec ted events. All events, however, were always counted with one of 24 accumulators ()Griin et al 



19951 ^. This way, the data can be corrected for incomplete data tr ansmission so that reliable impact 
rates can be determined for all satellite flybys ( Kriiger et al.l . 2000 ). 



Since its injection into an orbit about Jupiter, the Galileo spacecraft has been exposed to the harsh 
radiation environment of the Jovian magnetosphere with energetic particles of up to several MeV 
energies. The Galileo Jupiter mission was extended three times so that the spacecraft was exposed to 
a total radiation dosage five times higher than it was originally designed for. Especially high radiation 
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dosages were acquired during orbit insertion in December 1995 and during several perijove passages 
after mid-1999 when Galileo's perijove distance from Jupiter was about 6Rj. It was anticipated that 
these high radiation levels would cause severe damages to the spacecraft electronics and the scientific 
instruments. Although degradation of the dust instrument was recognised in the dust data, no failure 
has occurred so far. The degradation effects include - amoung others - a drop of the channeltron 
amplification, shifts of the measured instrument current, charge rise times and amplitudes which 
reduced the sensitivity for dust impacts and noise events (Kriiger et al., in prep.). The most important 
effect for our analysis here is shifts in the speed and mass calibration of the dust impacts. 



2.2 Impact direction 



The analysis of the dust measurements obtained at Ganymede showed that the impact direction of 
the particles could be used as one important parameter to identify ejecta particles belonging to a dust 
cloud surrounding this moon ( Kriiger et al.l . I2OOCI ) : in particular, the impact direction of the g rains 
could be used to separate cloud particles from Jovian dust stream particles ()Griin et al.l . Il998h . As 
rotation angle, Q, we define the viewing direction of the dust sensor at the time of particle impact. 
During one spin revolution of the spacecraft, the rotation angle scans through 360°. Rotation angles 
for the Galileo dust instrument, however, are reported opposite to that of the actual spacecraft rotation 
direction. This is done to easily compare Galileo results with the dust detector data taken on the 
Ulysses spacecraft, which, unlike Galileo, has the opposite spin direction. Zero degrees of rotation 
angle is taken when the dust sensor points close to the ecliptic north direction. At rotation angles 
of 90° and 270° the sensor axis lies nearly in the ecliptic plane (which is close to Jupiter's equatorial 
plane) . 

The dust instrument itself has a 140° wide field of view (FOV). Dust particles which arrive within 10° 
of the positive spin axis can be sensed at all rotation angles, while those that arrive at angles from 
10° to 130° from the positive spin axis can only be sensed over a limited range of rotation angles. A 
sketch of the detection geometry at close satellite flybys is shown in Fig. ^ 



2.3 Impact velocity 



Calibrated impact velocities are derived from the rise times of the impact charge signals by an em- 
pirically derived algorithm ( Griin et al.l . 1995 ). The analysis of the dust impacts detected close to 
Ganymede showed that their average impact veloci ty onto the detecto r target was close to the en- 
counter velocity of Galileo with this moon (8 kms~^) ( Kriiger et al.l . l2000l ). It implied that the particles 
truly originated from Ganymede and that they belonged to a steady-state dust cloud surrounding this 
satellite. This good agreement of the measured mean impact velocity with the expected velocity also 
showed that the calibration of the dust instrument is reliable in this velocity range. 



Two statistical subsets of particles could be separated: nearly all of the Ganymede particles had 
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calibrated velocities below 10 kms~^, whereas most of the stream particles had higher velocities. The 
calibrated impact velocity has been used as a parameter to separate both populations of dust at 
Galileo's G8 Ganymede flyby when the Jovian stream particles approached the dust detector from the 
direction towards Ganymede and particles belonging to Ganymede's steady-state dust cloud could not 
be identified by their impact d irection alone. The true impact velocities of stream particles exceeded 



200 kms"^ (Zooketal 



IDc 



119961 ^ and were much faster than the velocity range of the dust instrument 
calibrated in the laboratory (70 kms~^). Thus, the velocities for the stream particles derived from 
the instrument calibration significantly underestimate the true particle velocities. 



3 Data analysis 

3.1 Europa 

3.1.1 Impact direction 

The antenna of Galileo usually pointed towards Earth for data transmission. This fixed the spin axis 
of the spacecraft so that the detector basically scanned the anti-Earth hemisphere. In addition, due to 
the orbital motion of Jupiter about the Sun, the geometry for dust detection with the dust instrument 
gradually changed with time, leading to the non-detectability of dust particles in the close vicinity 
of the Galilean moons after mid-1999. For this reason, ejecta particles were measurable during ten 
close Galileo flybys at Europa out of 11 flybys in total. No data could be collected during two of 
these encounters due to spacecraft anomalies (safings) so that data sets from eight Europa flybys are 
available (E4, E6, Ell, E12, E14, E15, E17, E19; see also TablUl. The labels of the encounters are: 
the first letter of the satellite encountered by Galileo plus the number of Galileo's orbit about Jupiter. 

In Fig. 121 we show the impact direction (rotation angle) of the dust particles detected within about 
2 h around closest approach to Europa whose complete set of measured impact parameters has been 
transmitted to Earth. During most flybys, particle impacts with 180° < < 360° were concentrated 
towards Europa. This is most obvious during encounters E4, Ell, E12 and E19. Most of these impacts 
were detected at altitudes below 3 Re (Re is the Europa radius, see Table QJ). 

To analyse the impact direction of the dust grains onto the detector we assumed that the speed of 
dust relative to Europa in the vicinity of the moon is low compared to Galileo's flyby speed. Thus, 
the approach direction of the dust for an observer moving with the spacecraft is more or less parallel 
to the velocity vector of Europa relative to the spacecraft (the so-called ram direction). Since the 
orbital planes of Europa and Galileo about Jupiter coincide to within a few degrees, such particles 
approached the detector from a direction corresponding to ~ 270° rotation angle during all eight 
encounters. Rotation angles of about 90° are opposite to the direction towards Europa. With the 
sensor field of view of 140°, particles detected with rotation angles 180° < G < 360° are compatible 
with an origin from Europa itself. In the following we will call them Europa particles. This detection 
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geometry is very similar to the one at the majority of the Ganymede flybys ([Kriiger et alJ . 120001 ) 



Insert Figure [2] 



The direction from which the Jovian dust stream particles were observed varied during Galileo's path 
through the Jovian system: when Galileo approached the inner Jovian system, rotation angles around 
270° were observed. Between 1996 and early 1999 (the time span considered here) the rotation angle 
shifted to 90° shortly before Galileo's closest approach to Jupiter and the stream particles approached 
from this direction on the outbound portion of the spacecraft trajectory. Therefore, depending on 
when an individual satellite flyby occurred, stream particles approached the sensor from one or the 
other direction. In the cases of the Europa flybys considered here, the stream particles approached 
from rotation angles 0° < < 180° {i.e. opposite to that of the Europa particles) or stream particle 
impacts had already ceased because of the unfavourable detection geometry. It should be emphasized 
that during most of these eight Europa encounters more impacts were detected from the Europa 
direction than from the direction from which stream particles were to be expected. Only between 
zero and two impacts from the direction of the stream particles occurred during six encounters. Only 
encounters E6 and E17 showed the same number of impacts from 0° < < 180° (stream particles) as 
from the opposite direction (Europa particles). The statistics of particle detections is given in Table ITll 



Insert Table |lT| 



A total number of 64 Europa particles have been identified below 8 Re altitude during these eight 
encounters purely by their impact direction (Table For our further analysis we use a cut-off 
altitude of 8 Re because this is close to the extension of Europa's Hill sphere (Table P). We can 
minimize the potential contamination by particles belonging to other Jovian dust populations this 
way (see Sect. 13.1.31 for a discussion of their impact rates). For instance, a cut-off altitude of 10 Re 
would increase the number of grains classified as Europa particles by only five. 

For some flybys, the numbers of identified Europa particles are lower limits to the true numbers of 
detected grains because the complete set of parameters measured upon impact could be transmitted 
to Earth for only a fraction of all impacts (column 9 of Table ITlj) . At E14, E15 and E19, however, the 
complete set of parameters was transmitted for all impacts within 2 h around closest approach. 

In Table ^ (columns 7 and 8) we compare the number of Europa particles with the number of all 
events (dust plus noise) detected by the instrument below 8 Re and from a direction 180° < < 360°. 
This shows that the noise contribution to the total number of detected events in this altitude range 
is between and 60% with an average of 23% (17 out of 81 events are classified as noise). Although 
Europa orbits Jupiter in the region where high noise rates occurred, the total number of noise events 
in the data set is relatively small. A plot similar to Fig. |21 but with all detected events (column 8 of 
Table HI)) also shows dust concentrations at the Europa closest approaches so that the derived densities 
would also peak towards the satellite (Sect. IXT3)) . Hence, our conclusion that the dust impacts are 
concentrated towards the surface of Europa does not depend upon the applied noise removal algorithm. 



9 



We will come back to the noise problem in Sect. 13 .131 where we will determine the spatial distribution 
of dust surrounding Europa. For our analysis of the Europa dust cloud we will use class 3 and noise- 
removed class 2 data. It should be noted that there is no physical difference between dust impacts 
categorised into class 2 and class 3. 



3.1.2 Impact velocity 

In Fig. 121 we show the velocity distribution of th e Europa particles whose impact velocity has been 



determined with a velocity error factor VEF < 6 (IGriin et al.l . ll99a ). 53 particles fulfill this criterion. 



During all eight encounters of Galileo with Europa, the flyby velocity was close to 6 kms~^ (Table ITlj) . 
well above the detection threshold of the dust instrument for micrometer-sized grains at 2kms~^. We 
can separate two subsets of particles from their velocity distribution, similar to the measurements 
at Ganymede: the Jovian stream particles with velocities typically above 10 kms^^ and the slower 
Europa particles. The mean velocity of the 53 Europa particles in Fig. El is 5.5 it 3.5 kms~^ {^^)- 
Given a typical uncertainty for an individual velocity measurement of a factor of two, this value is in 
good agreement with the velocity of Galileo relative to Europa. 



Insert Figure |31 



The velocity measurements ~ like the measurements at Ganymede - are in agreement with dust 
particles belonging to a dust cloud of Europa. They confirm that the empirical velocity calibration 
of the dust instrument can be applied to the relatively slow ejecta particles, although the calibration 
is wrong for the much smaller and faster Jovian dust stream particles. This result also confirms that 
the calibrated impact velocities can be used to identify particles belonging to a dust cloud when such 
grains cannot be identified by their impact direction. This will be applied to the data obtained in 
the vicinity of Callisto in Sect. 13.2.21 It should be emphasized that the two velocity distributions in 
Fig. O overlap, leading to some ambiguity in the identification of individual grains. Hence, the velocity 
criterion can only be applied to a statistically large data set. 

In Fig. Owe have marked particles according to their calibrated velocities: those with impact speeds 
below lOkms^^ are shown as circles, faster grains as crosses. During seven Europa encounters of 
Galileo, the majority of particles with 180° < < 360° had impact speeds below 10 kms^^, consistent 
with a particle origin from Europa. Thus, for these seven flybys the identification of cloud particles 
with the velocity criterion is in agreement with the identificaton by the impact direction alone. Only 
the E19 encounter had 50% of particles (5 out of 10) with higher calibrated impact speeds so that 
the majority of grains would be classified as stream particles and be rejected if the impact direction 
were not applicable as the main identification parameter. Most of the particles detected at E19 with 
impact speeds above lOkms"^ (4 out of 6 particles with 0° < < 360°), however, have a velocity 
error factor VEF > 6, which makes their speed calibration very uncertain, anyway. To summarize, 
the identification of Europa particles from their impact direction and impact speed is quite reliable 
for all eight Europa encounters. 
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3.1.3 Impact rate and number density 



With 64 complete data sets of particles detected during eight Europa encounters we can calculate the 
dust impact rate in the close vicinity of this moon (Fig. ^ . We have defined distance bins equahy 
spaced on a logarithmic scale. Then we divided the number of particle impacts in each bin for which 
the complete set of measured impact parameters has been transmitted to Earth by the time Galileo 
has spent in that bin (dotted lines). Finally, to correct for incomplete data transmission (Sec. [2^, 
we have multiplied the impact rate bin by bin with the ratio between the number of counted particles 
and the number of particles for which the complete data set has been transmitted. These corrected 
impact rates are shown as solid lines. 



Insert Figure |4] 



Figure 13 shows that for those Europa encounters where the number of detections is sufficiently large 
(at least ten particles; E4, Ell, E12 and E19) the impact rate clearly increases towards Europa. This 
implies a concentration of dust particles at Europa. It wa s already obvious fro m Fig. |21 and confirms 
the earlier results of a dust concentration at Ganymede ((Kriiger et alJ . [20001). On the other hand, 
at the remaining four encounters, the number of detections is so low that no statistically meaningful 
radial profile can be derived although the data are also compatible with particle concentrations at 
Europa. It should be noted that the correction for incomplete transmission is small in all bins and 
does not significantly affect the slopes of the power law fits. The slopes derived with correction for 
incomplete transmission are in the range —1 .4 to —2.7 (Table [TTI) . This is somewhat fiatter than the 



most reliable slopes obtained for Ganymede (jKriiger et al 




20001 ^. 



An important question arises: is the slope of the dust distribution at Europa truly flatter than that at 
Ganymede or is it an artefact caused by the measurement process? Two effects may cause a flattening 
of the slope: 1) a background of dust particles in jovicentric space; and 2) incomplete removal of noise 
events by our noise removal technique. Both would lead to an artificial flattening of the derived dust 
impact rate profile. 

To analyze the first hypothesis - a background dust population - we consider the complete data 
set of Galileo dust measurements in the Jovian system: the dust instrument has detected a numbe r 
of micrometer-sized particles mostly in the region between the Galilean moons (|Grtin et al.l . |l99a). 
At least two populations of grains can be distin guished: part i cles on b ound prograde orbits about 



Jupiter and a popu lation on retrograde orbits l|Colwell et al 



Krivov et al 



1998 



Thiessenhusen et al 



200C 



2002a|). Depending on the detection geometry of the dust instrument during a specific 

3 oth p opulations taken together were up 



orbit of Galileo about Jupiter, impact ra tes of particles from 



to six per day in the region of Europa ((Thiessenhusen et al 



2000l . their Fig. 10). Considering that 



one Galileo passage through Europa's Hill sphere (Table QJ) lasted about 75 min, only one dust impact 
from these populations has to be expected every third Europa flyby. We therefore conclude that a 
contamination by particles on jovicentric orbits is negligible for our analysis of the Europa dust cloud. 



11 



A potential contamination caused by such grains would be even smaller at Ganymede and Callisto 
because their number density decreases further away from Jupiter. 

In order to check the second possibility - incomplete removal of noise events from the data set - we 
use an alternative approach to calculate noise-free dust impact rates. The noise removal algorithm 
applied so far determines whether each individual class 2 event is most likely a true dust impact or 
a noise event. Instead, we calculate the average noise rate measured with the dust instrument and 
subtract it from the total counted rate to obtain the dust impact rate. We have first calculated the 
total event rate of dust impacts plus noise with the same technique as before, i.e. from the complete 
data sets, without applying our noise-removal scheme. This gives somewhat flatter power law slopes 
than those derived for the noise-removed data set (between —1 and —2). We have then calculated the 
fraction of noise events in the class 2 accumulator data with our noise identification scheme during 
a one-day interval centered around each Europa flyby and calculated the rate of noise events in the 
counter data by multiplying the total counted rate with the fraction of noise events. Typical noise 
rates are between 0.04 and 0.1 per minute. We have then subtracted this noise rate from the total 
event rate obtained from the entire data sets. The resulting radial density profiles have power law 
slopes between —2 and —4. The theoretically expected value is about —2.5, or more precisely it is 
steeper than —2.5 for r ^ 5i?sat (bound gr ains dominate) and it becomes flatter, between —2.5 and 



2 (escaping grains dominate), farther out (jKrivov et al .1 . 120031 ). Thus, we conclude that the observed 



radial density profiles which are somewhat flatter than the expected values may be due to incomplete 
noise removal. 

We do not investigate variations of the slopes between individual encounters because of the large 
statistical uncertainties and the potential unrecognised noise contamination of the data. Spatial 
variations with respect to the flyby position relative to the satellite will be addressed in a future 
investigation. 



Insert Figure \S\ 



With the impact rate profiles derived for the individual Europa encounters we can now calculate the 
spatial density of dust in the environment of this satellite. We first divide bin by bin the impact rate 
by the spin-averaged detector area to obtain fiuxes (m~^s~^). Then we divide these fluxes by the 
mean impact velocity (spacecraft velocity relative to the moon) for a given flyby. This results in mean 
number densities (m~^) in the various distance bins. Note that the slope of the number density is the 
same as that of the impact rate, because both the spin-averaged detector area and the mean impact 
velocity are assumed to be constant (independent of distance) for any individual flyby. The result is 
shown in Fig. |5J The number densities show a clear increase towards Europa and the average slope is 
—2.02 ± 0.63. It is remarkable that the variation in the derived number densities from encounter to 
encounter is relatively small. Since the closest approaches of Galileo at Europa occurred at different 
longitudes and latitudes of Europa, it indicates that the dust distribution around this moon does not 
show strong spatial or temporal variations. 
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We have also checked alternatives for the origin of the dust impa cts detected in the close vicinity of 



the Galilean satellites other than the impact-ejection mechanism (jKriiger et all 1200(1 ) . Gravitational 
capture of the grains by the satellites can be dynamically ruled out. Electromagnetic interactions seem 
to be too weak, in particular at Europa and Callisto which do not have their own magnetic fields. 
Although the geysers on lo are the most likely source for the Jovian dust stream particles, no geyser 
activity has been observed on the other Galilean moons. The most plausible explanation for the origin 
of the dust grains was continuous bombardment of the satellites by interplanetary micrometeoroids. 

In this work we assume spherical symmetry of the clouds which is supported by our measurements: 
passages of Galileo at different latitudes and longitudes of the moons did not reveal strong asymmetries 
in the dust distribution, except, possibly, at Callisto (see Sect. 13.23)) . This implies a spherical structure 
of the dust distribution surrounding the satellites. Of course, it does n ot rule out the existen ce of 



asymmetries in the dust density which will be investigated in the future ([Sremcevic et al.l . 120031 ) . 



We now look at the number density profile expected from theory. iKrivov et al.l (|2003f ) developed a 
model of a spherically symmetric, stationary dust cloud around a satellite, maintained by impacts of 
interplanetary micrometeoroids. To the first approximation, the number density of dust grains ejected 
into ballistic orbits, which dominate the cloud at distances of several satellite radii, is 

rihoundix) oc x"^/^, (2) 

where x = r/Rs^t is the distance measured in satellite radii. The contribution of escaping grains into 
the cloud is somewhat shallower: 

nunbound(a;) oc (3) 

which slightly flattens the overall radial profile at larger distances from the moon, closer to its Hill's 
sphere. 



Krivov et al 



( 20031 ) have also constructed an algorithm to calculate the proportionality factors in 
Eqs. (HI)-®- The algorithm implies a chain of estimates: for the mass flux and typical speed of 
projectiles, for gravitational focussing of impactors by Jupiter, for the ejecta yield, ejecta mass and 
velocity distributions, etc. The values of the model parameters (both assumed and derived) for 
Europa are given in Table ITm Other parameters that have the same values for all Galilean satellites, 
are: slope of cumulative ejecta mass distribution a = 0.83, maximum mass of an ejected fragment 
Afmax = 10~^g, opening angle of cone into which particles are ejected ipo = 90°. For a detailed 
description of the parameters, the reader is referred to the original paper. We note that the model 
calculates number densities of particles with masses above the detection threshold of the Galileo dust 
detector. As the threshold is speed-dependent, the number densities are computed separately for each 
Galileo flyby, and the results for specific flybys generally differ even for the same moon and the same 
distance. A strong advantage of this approach is that it enables direct comparison of the number 
densities predicted by the model with those derived from the measurements. 

Using this model for Europa, we obtained theoretical curves superimposed on the data points in 
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Fig. [51 A comparison between the number densities derived from the Gahleo measurements and those 
computed with the model will be given in Sect. |3 

Insert Table HITI 



3.1.4 Mass distribution 



The charge released by an impact of a dust particle onto the detector target depends on the mass and 
the velocity of the grain (Equ.^). In particular, to calculate the particle mass one has to know its 
impact velocity. The calibration of velocity and mass from the measured charge rise times and charge 
amplitudes is usually performed based on laboratory measurements obtained at a dust accelerator. 

In Fig. ini we show the mass distribution of the particles from all eight Europa flybys for which the 
velocity could be reliably determined (VEF < 6; 53 particles in total). In the upper panel the 
complete instrument calibration has been used to obtain particle speed and mass. With this method 
the uncertainty of the impact velocity is typically a factor of 2 and that of the mass is a factor of 10. 



Insert Figure El 



The dust detector has a velocity-dependent detection threshold (|Grun et al.l . ll99fJ ). The threshold for 
particles approaching with 6kms~^ is shown as a dashed line. The mass distribution is incomplete 
around this value. 

The mass distribution is also affected by the low data transmission capability of Galileo and the data 
storage scheme in the instrument memory. As a result, nearly all data sets lost are in the lowest 
amplitude range ARl which — for particle velocities of about 6kms~^ — corresponds to the mass 
range below ~ 3 X 10"^^ kg. If we assume that the lost particles are equally distributed over the 
mass bins below this value, the maximum of the mass distribution is artificially too low by less than 
a factor of 1.2. Thus, incomplete data transmission does not significantly affect the mass distribution 
for Europa particles. 

If the individual impact velocities of dust particles were known with a higher accuracy than the typical 
factor of 2 uncertainty from the instrument calibration, the uncertainty in the mass determination 
could be improved. The measured mean impact velocities of Europa particles are close to the velocity 
of Galileo relative to Europa during the individual encounters (Table ^) . We therefore assume the 
latter ones as the particles' impact velocities and show the recalculated particle mass in the lower 
panel of Fig. (21 The width of the mass distribution is significantly smaller than that derived from 
the calibrated impact velocities. This me thod has also been successfully applied to calculate the size 
distribution of Ganymede eiecta particles ( Kriiger et al. , bflOC ) and interstellar dust particles measured 



with Galileo and Ulysses ([Landgraf et al.l . 120001 ). 

The mean mass of the Europa particles is 9.0 x 10^^^ kg. Assuming spherical particles with a density 
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of 1 g cm ^ - the density of water ice - this corresponds to a particle radius of ~ 1 ^m. 

Degradation of the dust instrument caused by the high radiation dosages in the Jovian magnetosphere 
lead to shifts of the calibrated masses and impact speeds of the dust particles. All data collected 
after mid-1997 are affected by this shift, the later in the mission the data were collected, the stronger 
the shift. For our Europa measurements this means that masses are too low by a factor of about 1.5 
beginning with the Ell encounter. We have corrected the calibrated masses for these data sets and 
constructed a corrected mass distribution (solid histogram in the bottom panel of Fig. Ej). Since we 
have taken the speed of Galileo relative to Europa as the impact speed of the particles, we need to 
correct the masses only. The shift in the velocity calibration caused by the instrument degradation does 
not affect this mass distribution. The resulting mean mass of the Europa particles is 1.3 x 10^^^ kg. It 
should be emphasized that the shift in the velocity calibration does not affect the identification of the 
Europa particles because the particles at the Europa encounters have been identified by their impact 
direction alone, without using the impact speed as an additional criterion. 



3.2 Callisto 

3.2.1 Impact direction 



During Galileo's prime mission about Jupiter in 1996 and 1997 the spacecraft had three close flybys 
at Callisto (C3, C9, CIO). The spacecraft orientation during these encounters allowed the detection 
of ejecta cloud particles close to Callisto. During all Callisto flybys after mid-1999 the spacecraft 
orientation prevented the detection of ejecta cloud particles so that the measurements at Callisto are 
restricted to these three encounters. Unfortunately, all three of them occurred on the portion of the 
Galileo trajectory inbound to Jupiter where Jupiter stream particles and potential ejecta particles 
from Callisto approached the dust sensor from the same direction (180° < < 360°; Sect. I3.1.1|l . 
This is shown in Fig. [7| almost the entire number of dust impacts measured close to Callisto occurred 
from this direction. Thus, Callisto particles could not be uniquely identified by their impact direction 
alone, and we had to use the impact velocity as an additional criterion to identify them. 



Insert Figure [71 



The analysis of the ejecta particles detected at Europa (Sect. I3.1.2() and Ganymede (|Kriiger et al 



2000|) showed that - on average - stream particles and ejecta cloud particles occupy different regimes in 
calibrated impact speed: cloud particles have typical speeds below 10 kms^^ which are on average very 
close to the encounter velocity of Galileo with the satellite, whereas stream particles have significantly 
higher calibrated speeds. This has been successfully applied to identify Ganymede cloud particles from 
Galileo's G8 encounter at Ganymede. We apply the same velocity criterion here to separate Callisto 
particles from the dust streams. The numbers of Callisto particles identified this way are listed in 
Table Ull for each orbit. The total number of Callisto particles from all three encounters is 35. 
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Only particles detected at altitude below 6 Rc are considered for further analysis in order to minimize 
the contamination by stream particles. The analysis of the Europa and Ganymede data showed that 
a few stream particles have calibrated velocities below 10 kms^^ and would erroneously be classified 
as Callisto particles (Fig. EJ. Most clo ud particle irnpacts at Europa and Ganymede occurred below 
an altitude of about 6Rsat (Fig. |2 and ( Kriieer et alJ . l200(t )) so that the inclusion of particles detected 
further away would increase the probability that the particles are actually stream particles rather than 
dust cloud particles. An apparent concentration of stream particle impacts within 3Rc altitude at C3 
and within 5Rc at CIO (Fig. [TJ is due to a higher data transmission rate of Galileo in these periods 
( Kriitrer et all . boOlI 'l. 

Callisto orbits Jupiter outside the region within 15 Rj where the high noise rates occurred. Thus, a 
potential noise contamination of the Callisto data is expected to be very low. For our analysis of the 
Callisto dust cloud we will use class 3 and class 2 data without noise removal. 



3.2.2 Impact velocity 

A total number of 35 Callisto particles have been identified by their calibrated impact velocity below 
10 kms~^ and below an altitude of 6Rc. Their mean impact velocity is 6.4 it 2.1 kms~^ (l*^)- This 
value is artificially too low because the velocity distribution is cut off at 10 kms~^. The average flyby 
speed of Galileo at Callisto was 8.1 kms~^ and both speeds agree within 1 a. The velocity distribution 
of the Callisto particles is shown in Fig. |HJ 

Insert Figure |H1 

Degradation of the dust instrument electronics does not significantly affect the impact velocities of the 
particles because most the Callisto encounters occurred relatively early during the Galileo mission. In 
particular, the identification of Callisto particles via their impact speed is not affected. 

3.2.3 Impact rate and number density 

With Callisto particles from all three Callisto encounters identified by their impact speed and impact 
direction, we can construct the radial profile of the dust impact rate in the same way as we did for the 
Europa flybys before. This is done in Fig. El For the C3 and CIO encounters the impact rate increased 
towards Callisto. At C9 the number of detected particles was only three. The derived radial profile 
— although being very uncertain — is compatible with an increase towards Ganymede. We conclude 
that the dust impact rates detected at Callisto are compatible with a dust concentration surrounding 
this moon. 

Insert Figure El 

In Fig. ^1 we show the number densities for Callisto derived from the radial profiles of the impact 
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rate. The data from the individual flybys show an increase towards the surface of the moon. The 
number densities derived from the three flybys show a large variation from flyby to flyby, much larger 
than those for Europa. In particular, the number densities derived from the C9 data are very low. 
Possible reasons for this variation will be discussed in Sect. IH 

In the same figure, the su perimposed curves show the number density profile calculated with the 



model (Krivov et al 



2003) and parameters listed in Table ITllI A comparison between the data and 



the model will be given in Sect. |3 



Insert Figure 1101 



3.2.4 Mass distribution 

The mass distribution for the 35 Callisto particles with VEF < 6 is shown in Fig. 1111 As for Europa, 
we show the mass distribution with the calibrated impact velocities (upper panel) and that obtained 
by applying the velocity of Galileo relative to Callisto (bottom panel). Here the detection threshold is 
8kms~^, and the mass distribution is incomplete around the threshold. Again, the mass distribution 
is incomplete in the bins below about 10~^^ kg due to incomplete data transmission and the maximum 
of the mass distribution may be artificially too low by up to a factor of 1.3. 

Insert Figure 1111 

The mean calibrated mass of the Callisto particles is 3.7 x 10~^^kg, which is an order of magnitude 
lower than the value derived for the Europa cloud particles. Again, assuming spherical particles 
with a density of lgcm~^, this corresponds to a particle radius of 0.5 /xm. We have also corrected 
the masses of the particles for instrument degradation and constructed a corrected mass distribution 
(solid histogram in the bottom panel of Fig. Ill|). The corrected mean mass of the Callisto particles is 

5.2 X 10-1^ kg. 

3.3 lo 

3.3.1 Impact direction 

Galileo had a total of seven flybys at lo but only the initial one in December 1995 (10; the orbit 
notation is I "zero") had a favourable detection geometry. During the other flybys at this satellite 
which occurred after mid-1999 the sensor orientation prevented the detection of ejecta dust particles. 
Figure [T^ shows the sensor orientation at particle impact at the 10 encounter. The bottom panel shows 
class 3 and noise-removed class 2 data (only four impacts). 

The noise identification criteria applied to the Europa data have been developed for the spatial region 
outside 10 Rj. The 10 data, however, have been collected closer to Jupiter at 6Rj where the noise 
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characteristics may have been different (jKriiger et all Il999d ). We therefore show the full data set of 
classes 1 to 3 in the top panel of Fig. 1121 The noise fraction in class 2 derived with the algorithm 
for secondary ejecta grains is about 80%. The complete class 2 data set also shows a concentration 
of grains towards lo. Therefore, the noise rejection algorithm may be too restrictive, thus rejecting 
too many events. Class 1 events, which are usually classified as noise in the Jovian environment, 
show an interesting behaviour: the impacts cluster at rotation angles < © < 180°. This direction is 
compatible with the approach direction of plasma particles in the lo torus. 

A few days before the fl yby at lo, the channeltron high voltage w as decreased and the charge detection 



thresholds were raised ( Griin et al 



1996 



Kriiger et al 



Il999bh 



to reduce the instrument sensitivity 
in the high radiation environment of the inner Jovian magnetosphere. This reduced the instrument 
sensitivity for class 3 impacts. Those impacts, however, that did not generate enough charge to become 
class 3 events should have shown up as class 2 impacts. Unfortunately, class 2 is contaminated with 
noise so that the identification of these particles is ambiguous. Fortunately, only a small fraction of 
the data sets of particles was lost due to incomplete data transmission (Table Hi]) . 



Insert Figure 1121 



Figure El shows four particles within an altitude of about 6 Rj (the Hill sphere of lo) . Two of these 
particles approached the detector from a direction ~ 270°. The other two particles were detected 
when the dust detector pointed ~ 90° away from this direction. A check of the approach direction of 
potential lo particles revealed that all four particles are compatible with an lo origin: the approach 
direction of lo particles as seen from Galileo was so close to the spin axis of th e spacecraft (directio n 
opposite to antenna axis) that they were detectable at almost all rotation angles (|Kruger et al.l . ll999d ). 
We therefore consider all four particles as probable lo particles. Only one of these impacts occurred 
in ion amplitude range ARl, two were d etected in AR2 and one in AR3. lo dust stream particles were 
detected in ARl only ( Griin et al. . 1998l l so that the identified particles were unlikely stream particles. 
In addition, the impact rate of strea m particles was reduced in the inner Jovian system because of 
reduced impact speeds (|GraDsl . l200ll ^. 



3.3.2 Impact velocity 



The total number of identified potential ejecta particles from lo is only four. With this very low 
number of detections no reasonable velocity distribution can be constructed. We can, however, still 
check whether the average particle velocity is compatible with the hypothsis that the particles are 
basically at rest with respect to lo. The velocity of Galileo relative to lo at the 10 encounter was 
15kms^^. The averaged particle speed derived for the four particles is 10.3 it 8.4 km s^^. Although 
one has to keep in mind that the statistical uncertainty of this value is very large, it is compatible 
with particles being bound to lo. 
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3.3.3 Impact rate and number density 



The dust impact rate at lo derived from the four identified lo particles is shown in Fig. 1131 It shows 
a sUght concentration towards lo. One has to keep in mind, however, that the radial profile of the 
impact rate is very uncertain because of the small number of dust detections in the close vicinity of 
this moon. 



Insert Figure 1131 



Insert Figure 1141 



The number density of dust in the close vicinity of lo derived from the impact rate profiles is shown in 
Fig. 1141 For comparison we show the data points for the noise-removed data set (solid lines) and for the 
complete class 2 and class 3 data set. The curves in the same figure depict the number density profile 



calculated with the model (Krivovetal 



20031 ) and parameters listed in Table III II A comparison 



between the data and the model will be given in Sect. 0] 



3.3.4 Mass distribution 



No attempt has been made to construct a mass distribution because of the small number of detections 
at this moon. The mean mass of the particles taking the speed of Galileo relative to lo as the impact 
speed is 8.5 x 10~^^ kg. Note that this value is not affected by the instrument degradation because 
this lo flyby occurred at the beginning of Galileo's Jupiter mission. 



4 Comparison of the dust clouds surrounding the GaUlean moons 



In the previous section we have analyzed the dust impacts detected in the circumsatellite dust clouds 
individually for lo, Europa and Callisto . The dust cloud of Ganymede has been investigated in an 
earlier publication (|Kriiger et al.l . |2000|). For each of the four moons we have identified impacts of 
probable ejecta cloud particles in the Galileo dust data set, determined their impact speeds and mass 
distributions and have derived impact rate and number density profiles. We now take the data sets 
for all four moons together to compare the properties of their surrounding dust clouds. 



4.1 Mass distribution 



The mass distribution of the grains allows a simple check for the compatibility of the data with 
the hypothesis of the impact origin of the detected p articles. We t ook t he mass distributions for 
Europa, Callisto and Ganymede (Fig. I^J Fig. ^2 ^-nd (jKriiger et al.l . [20001)) and show linear fits to 
the cumulative distributions in Fig. 1151 The slopes of the cumulative mass distributions assuming 
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Galileo's velocity relative to each moon as the particle impact speed are given in Table HVl Only one 
data point is shown for lo because of the scarcity of the data. For the three other moons the slopes are 
in the range between 0.58 and 0.86, w hich is in good agreement w ith the typical slopes one expects for 
impact ejecta (0.5 ^ a ^ 1.0; see, e.Q.. lKoschnv and Griinl ()200ll )). It should be emphasized that even 
though the statistical uncertainties in the data sets are relatively large because of the small number of 
detections, the slopes derived for the three moons do not differ very much. The slopes derived for the 
mass distributions taking the calibrated impact speeds are in the range 0.5 - 0.6 and are thus even 
closer together, although they are flatter than those derived with the spacecraft speed relative to the 



moon. 



It should be noted that the mass distributions of dust around Ganymede and Callisto agree very well 
(0.82 vs. 0.86) whereas the one for Europa is somewhat flatter (0.58). This might reflect differences in 
surface properties of the satellites: for instance, flatter mass d istributions are typical of looser targets 



than of consolidated ones (see, e.g.. iKoschnv and Griinl l200ll . and references therein) 



Insert Figure 1151 



Insert Table ITVl 



4.2 Number densities 



The number densities derived for all four Galilean moons are shown in Fig. 1161 Straight lines are 
least squares fits to the data for each moon. We list the slopes of these curves, which are averages 
of slopes for individual flybys at each moon, in Table HVI (col. 8). The average slopes for each moon 
are between —1.6 for Callisto and —2.8 for Ganymede, with Europa being in between (—2.0). The 
very uncertain slope for the lo data (—2.0) is close to the Europa value. The Ganymede data show 
the steepest slope but also the largest uncertainty. This is mainly caused by the incomplete data 
transmission which mostly affected the G2 and G8 measurements. The incompleteness affects the 
error via oc (A^ it \/]V) x correction, therefore giving larger limits compared to the 100% transmission 
case. 



Altogether, the slopes ar e roughly con s istent with the one predicted in the framework of the spherically- 
symmetric cloud model (jKrivov et al.l . l200,'j ) . between —2.0 and —2.5 [see Eqs. ©-((SI)]. Only the slope 
derived for Callisto is somewhat flatter. A study of possible asymmetry effects in the clouds has shown 
that this "reference" slope may be substantially flatter or steeper , depending on the fiyb v geometry 
and the position of the satellite in its orbit at the time of flyby ((Sremcevic et al.l . 120031 ) . This may 
account for a tangible scatter in the slopes that we derived from the data. 

The measured absolute number densities of all four clouds at a given distance (measured in satellite 
radii) are similar — see Fig. 1161 It is important to compare this result with the theory. The absolute 
number density of a dust cloud at a given distance from a satellite should depend in a non-trivial way 
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on a number of factors: mass of the moon, its planetocentric distance (through a distance-dependent 
gravitational focussing of the impactor flux), as well as the satellite surface properties. All these 
depende nces are taken into account in the model which was confined, however, to a solid ice-silicate 
surface (|Krivov et all 1200,'^ ) . The number densities of the clouds around Europa, Ganymede and 
Callisto, computed with the model, turned out to be within a factor of several from each other (most 
notably, larger yields Y for closer-in satellites are compensated by lower ejecta speeds no, because of 
an energy conservation requirement used in the model). The number densities of these three clouds 
are consistent with the densities derived from the data. Not so for the lo cloud: the same model 



(Krivovetal 



2003) predicts a much higher dust number density (mostly because of the larger flyby 



speed at lo, resulting in a possibility of detecting much smaller dust grains compared to the other 
Galileans) . The number density in the lo cloud predicted by the model is at least an order of magnitude 
higher than observed. This could be due to the scarcity of the lo data (4 individual dust impacts 
only). Alternatively, this may be a real effect, caused by different surface properties of lo (volcanic 
deposits and condensed gases like SO2 frost) compared to the other three Galilean moons (presumably 
"solid" ice with some contamination by non- volatile materials). Such a view seems to be indirectly 
supported by our preliminary analysis of the dust en vironment between th e orbits of lo and Europa 
(work in progress): the density of the "Galilean ring" (jKrivov et al.l . l2002al ^ does not seem to increase 
from the Europa orbit towards Jupiter, which might be compatible with lo being a weaker source of 
ejecta than it would be, if it were similar in surface properties to the other Galilean satellites. 



Insert Figure 1161 



A comparison of data from different flybys at the same moon shows that for Europa the derived dust 
densities do not vary significantly from flyby to flyby. This indicates little or no temporal variation 
and/or dust density variation between the leading and the trailing side of this moon with respect to 
the field of impactors. On the other hand, the three Callisto flybys showed a significant variation by 
more than an order of magnitude between the C9 and the CIO flybys. For Ganymede, the G2 and G8 
flybys showed somewhat larger number densities than the Gl and G7 flybys. The former, however, 
have the largest uncertainties because of incomplete data transmission and particle identification via 
the impact speed criterion. At present, it is not clear whether the differences between the data from 
different flybys of Callisto and Ganymede c an be attributed to asymmetries in the circumsatellite dust 

I200J. A comparison between the data and theory is hampered 



clouds modelled in 



Sremcevic et al 



by poorly known directionality of impactors in the vicinity of Jupiter. This issue will be the subject 
of a future investigation. 

4.3 Grain velocities in the clouds 



The theoretical models of dust clouds predict that, on average, the co nstituent part i cles sh ould have 
substantial velocities with respect to their parent satellites. Based on iKrivov et al.l (|2003l ) formulas. 



about 40-50% of the grains at distances from 2 to 8 i?sat in the Ganymede and Callisto clouds should 
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be faster that 2kms ^. For Europa the fraction is 30-40%. 

The question that we address now is: can we find indications in the data that some grains have 
appreciable velocities relative to the respective moon? The most natural way would be to look at 
possible deviations of the actual impact speeds from the mean value, equal to the spacecraft speed 
with respect to the satellite. Unfortunately, this is not possible: as we have seen (Figs. OlandlHl), 
the instrument calibration is by far not accurate enough to do that. Another possibility would be 
to look at the impact directions, i.e. at the rotation angles of impacts. For most of the flybys the 
cloud particles, if they were at rest with respect to the moon, could only be detected in the rotation 
angle range = 270° it x, where the semi-width x of the detectability range is a (known) function 
of the FOV opening angle and the angle between the Galileo ram direction and its spin axis, /?. The 
semi-width x is 180° for (3 < 10° and decreases to 67° at /3 = 90°. Should a particle, which we identify 
as a cloud particle, have had a @ value somewhat outside the range Q = 270° it x, this would be an 
indication that the particle had an appreciable velocity, so that the impact velocity deviated markedly 
from the anti-ram direction. 

One should not expect the number of such events to be high: even very fast grains can only show up 
in this test if the direction of their velocity vector is appropriate, and only for some flyby geometries. 
We have checked all cloud particles and found four individual impacts of this kind: one in G7, one in 
Ell and two in E12. 



4.4 Mass budget 



As was the case for Ganymede earlier ((Kriiger et al.l . 120001 ). we give some general estimates con- 



cerning the mass bu dget of the dust clouds of the other Galilean moons as derived from the model 



( Krivov et al 



2003). The results, including new estimates for Ganymede, are collected in Table IVl 
Note that these are only crude estimates which are uncertain by at least one order of magnitude, per- 
haps even more. The expected steady-state masses of the clouds range from about 10 tons for Callisto 
to about 200 tons for lo; the lo cloud may, however, be lighter — see discussion in Sect. 14.21 Inter- 
estingly, the mass injection rate of the material into circumjovian space is similar for all four moons 
and is, in turn, comparable to the mass flux of impactors onto respective satellites, ~ 100 gs~^. This 
means that each satellite "redirects" nearly as much dust into the circumplanetary space as it receives 
from the interplanetary one. Of course, the mass/size and velocity distributions of the "incoming" 
and "outgoing" matter are generally quite different. 



Insert Table |Vl 
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5 Conclusions 



We have examined the dust impacts registered by the Gahleo dust detector in the immediate vicinity 
of lo, Europa and Calhsto during a total of 12 flybys at these Jovian moons. By analyzing impact 
directions and velocities and the ma ss distribution, as w ell as spatial locations of the dust impacts in 
comparison with model predictions ()Krivov et all 1200,'! ) , we have shown that the particles originated 
from the moo ns. Our analysis tech nique was similar to an earlier investigation of dust data collected 
at Ganymede (jKriiger et al.l . 1200(1 ) . The dust impacts recorded at all four moons are compatible with 
impact debris produced by hypervelocity impacts onto the surfaces of these moons. For the icy moons 
Europa, Ganymede and Callisto, the mass distributions of the detected grains, as well as the spatial 
dust densities derived from the measurements are in fairly good agreement with the predictions from 
the model of hypervelocity impacts of interplanetary dust particles (IDPs), assuming contemporary 
models of IDP flux at a heliocentric distance of Jupiter and a low-temperature ice-silicate target. For 
lo, the number of dust detections is too small to derive a reliable mass distribution. The number 
density obtained for this moon is more than an order of magnitude lower than predicted by the model 
which assumes a solid ice surface. The lack of detections may be due to a softer, fluffier surface of lo 
compared to the three icy Galilean moons. lo's surface is (at least partially) covered with volcanic 
deposits. 

This work continues the analysis of the dust clouds surrounding the Galilean moons and confirms 
the previous scenario of ejecta dust clouds generated by hypervelocity impacts of micrometeoroids. 
Up to now, this had only been tested at Ganymede. Our theoretical description is based on the 
physical conditions in the Jovian system as well as available laboratory data of hypervelocity impacts. 
We have neglected any spherical asymmetries of the dust clouds surrounding the moons. To a first 
approximation, this is supported by the data, especially the measurements taken at Europa. On the 
other hand, the Callisto and possibly Ganymede measurements show a variation between different 
flybys which might be indicative of a leading-trailing cloud as ymmetry caused by the motion of the 
moons through the field of impactors, as predicted by theory (|Sremcevic et al.l . l2003l ). This will be 
addressed in a future analysis. 

Most of the dust ejected from the surface is launched into bound orbits and falls back to the moon. 
These short-lived, but continuously replenished grains form the ejecta dust clouds. A tiny fraction of 
impact debris is ejected at speeds sufficient to escape from the moon entirely. The ejected mass is 
comparable with the incoming flux of IDP impactors. The escaping grains go into orbit about Jupiter 



and most of them will eventually be swept up by one of the Ga 



forms a tenuous dust ring surrounding the planet l|Krivov et al 



ilean sa tellites. A tiny fraction of them 



2002aD . This ring is by far too tenuous 



to be detected optically. By the impact ejecta mechanism, moons turn out to be efficient sources for 
dusty planetary rings. In particular, Jupiter's gossamer ring and Saturn's E ring are thought to be 
maintained by ejecta particles from smaller moons which orbit their parent planets within the rings. 
In November 2002 Galileo traversed the gossamer ring and had a close flyby of Amalthea, one of the 
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small Jovian moons which orbits the planet within the ring region. The dust measurements collected 
during this passage may give new insights into the dynamics and feeding mechanism of this dusty ring 
and about the significance of small moons as sources of dust. 

All celestial bodies without gaseous atmospheres (asteroids and planetary satellites of all sizes) should 
be surrounded by an ejecta dust cloud. The dust particles in the cloud are composed of surface 
material from the parent body and, hence, carry information about the surface from which they have 
been kicked up. Our analysis of the Galileo in-situ dust data has shown that spacecraft measurements 
near celestial bodies — which act as sources of dust — can be used as a new diagnostic tool to analyze 
the surface properties of these bodies. This is of particular interest for the Cassini mission which 
will investigate the Saturnian system beginning in 2004. The Cassini dust instrument will be able to 
measure the chemical composition of particles in the dust clouds surrounding the Saturnian moons. 
This way, the surface properties of the source moons can be investigated remotely. Interestingly, the 
in-situ dust measurements turn into a remote sensing technique where the dust instrument is used like 
a telescope for surface investigation. 
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Figure 1: Galileo's trajectory and geometry of dust detection during the E4 Europa flyby. The Galileo 
spacecraft is sketched in an orientation it was in during the flyby (see text for details). The directions 
to Jupiter, Earth and Sun are shown. C/A indicates closest approach to Europa, FOV the field of 
view of the dust instrument. The orientation of the dust instrument shown corresponds to a rotation 
angle 6 = 270°. At 90° rotation angle it points in the opposite direction. 
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Figure 2: Sensor direction (rotation angle, 0) versus altitude of the Galileo spacecraft above the 
surface of Europa at the time of dust impact. Data are shown for all eight Europa encounters during 
which data were successfully collected (E4, E6, Ell, E12, E14, E15, E17, E19). The radius of Europa, 
Re is given in Table Q The altitude range shown corresponds to a time interval of 2 h. Each symbol 
indicates a dust particle impact, and the size of the symbol indicates the impact charge created by 
the particle (10~^^ C < Qi < 10~^ C). Circles show particles with impact speeds below lOkms"^ and 
crosses show particles with higher speeds. Galileo did not traverse the region between the vertical 
dashed lines. In E19 no data were collected later than 19 min after closest approach ( corresponding 



to 4. 5 Re altitude) because of a spacecraft anomaly. Noise events have been removed ()Kriiger et al 



2nnil . their Table 4). 
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Figure 3: Impact speeds derived from the instrument calibration for dust particles detected below 
8 Re altitude at eight Europa encounters. The solid line shows the distribution for Europa particles 
(180° < e < 360°) and the dotted hue that for stream particles (0° < 6 < 180°). The mean impact 
speed is 5.5 it 3.5 kms~^. Only particles with a velocity error factor VEF < 6 
been considered (53 Europa particles). 
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Figure 4: Impact rates of dust particles detected during eight Europa encounters. The dotted 
histogram bins show the impact rates derived from the number of particles in each bin for which their 
complete information has been transmitted to Earth. The solid histograms show the same rates, but 
corrected for incomplete data transmission. In cases where only a solid line is visible, no correction 
was necessary because the complete information of all particle impacts is available in that bin. The 
vertical dotted lines indicate the minimum altitude reached by Galileo at closest approach. Error bars 
denote the \/iV statistical uncertainty, with N being the number of particles for which the complete 
information has been transmitted. Dashed lines are power law fits to the corrected impact rate 
(Table Hlj) . In E17 data were collected with a low rate of one instrument readout every 7 min which 
leads to a larger uncertainty in the impact time. Therefore, the impact time of the dust particles was 
shifted by -3.5 min to compensate for the reduced time resolution. At E19, no data were collected 
beyond 4.5 Re after Europa closest approach and the impact rate shown was multiplied by two outside 
this region to empirically correct for the missing data. 
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Figure 5: The number density of dust as a function of distance f rom the center of Eu ropa, derived from 
the data (symbols with error bars) and predicted by the model ( Krivov et al. . 2003I . lines). Horizontal 
bars for the data symbols indicate distance bins which were used in processing the data (see text 
for details), whereas vertical ones reflect \/iV errors due to a limited number of impacts. Since dust 
impacts are measured above a certain speed-dependent detection threshold of the instrument, different 
flyby speeds imply different minimum mass of the detected grains and therefore should give different 
dust number densities. This has been taken into account in the model. Two model curves are shown: 
for the slowest (Ell) and the fastest flyby (E14). 
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Figure 6: Mass distribution of the particles detected during eight Europa flybys. The upper panel 
shows the distribution obtained by using the measured impact speeds derived from the instrument 
calibration. In the lower panel the speed of Galileo relative to Europa has been assumed as the impact 
velocity in order to calculate the particle mass. The dotted histogram is without aging correction, 
the data for the solid histogram has been corrected for aging of the dust detector electronics. The 
vertical dashed lines indicate the detection threshold for particles which approach the detector with 
the velocity of G alileo relative to E uropa (about 6kms~^). Only the 53 particles with a velocity error 
factor VEF < 6 (jGriin et al.1 . Il99,lh have been cori^dered. 
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Figure 7: Sensor direction (rotation angle, G) versus altitude of the Galileo spacecraft above the 
surface of Callisto at the time of dust impact. Data are shown for all three Callisto encounters at 
which the spacecraft orientation allowed the detection of cloud particles (C3, C9, CIO). The radius 
of Callisto, Rc, is given in Tabled The altitude range shown corresponds to a time interval of 2 h. 
Each symbol indicates a dust particle impact, and the size of the symbol indicates the impact charge 
created by the particle (10"^'' C < Qi < 10~^^ C). Circles show particles with impact speeds below 
lOkms"^ and crosses show particles with higher speeds. Galileo did not traverse the region between 
the vertical dashed lines. 
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Figure 8: Impact speeds derived from the instrument calibration for dust particles detected at all 
three Callisto encounters below 6Rc altitude with 180° < G < 360° and calibrated impact velocities 
below lOkms"^ f solid line) and a bove lOkms^^ (dotted line). Only particles with a velocity error 
factor VEF < 6 (|Griin et al.1 . Il99,ll ) have been considered (35 particles with lOkms ). 
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Figure 9: Impact rates of dust particles detected during three Callisto encounters. Same notation as 
in Fig. 01 Only particles with impact speed < lOkms"^ have been considered and no noise removal 
has been applied. 
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Figure 10: Same as Fig. El but for Callisto. The theoretical profiles for C9 and CIO are very close to 
the one which is plotted for C3. 
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Figure 11: Same as Fig. El but for three Callisto flybys. The detection th reshold is shown for 8 km s ^. 
Only the 35 particles with a velocity error factor VEF < 6 dOriin et al.1 . [l995) have been considered. 
The dotted histogram is without aging correction, the data for the solid histogram has been corrected 
for aging of the dust detector electronics. 
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Figure 12: Sensor direction (rotation angle, G) versus altitude of the Galileo spacecraft above the 
surface of lo at the time of dust impact for the 10 flyby. The radius of lo, Ri is given in Table d 
The altitude range shown corresponds to a time interval of about 45 min. Each symbol indicates an 
impact event (dust or noise), and the size of the symbol indicates the generated charge (2 x 10~^^ C < 
Qi < 10~^ C). Circles show class 2 and class 3, asterisks show class 1 in the lowest amplitude range, 
crosses denote class 1, amplitude range 2 and higher. Top panel: all data in classes 1 to 3, bottom 
panel: class 3 and noise-removed class 2 data. 
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Figure 13: Impact rates of dust particles detected during Galileo's 10 encounter of lo. Same notation 
as in Fig. |11 Solid histogram: class 3 and noise-removed class 2 data (derived from only four dust 
impacts), dashed histogram: class 3 and all class 2 data (without noise removal). 
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Figure 14: Same as Fig. [S] but for To. Class 3 and noise-removed class 2 data are shown with solid 
lines. All class 3 and class 2 data taken together are shown with dotted lines. 
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Figure 15 : Cumulative mass distributions from Fig. El and ^2 for Europa and Callisto, respectively, 
and from iKriiger et 

ID diooo) for Ganymede. For lo only one data point is shown because of the 
sparcity of the data. Straight lines are linear fits to the data. Vertical bars indicate the ^/N statistical 
uncertainty. 
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Figure 16: Number density of dust as a function of radial distance from the center of Ganymede 
(data from 4 flybys), Europa (8 flybys), Callisto (3 flybys) and lo (1 flyby). The altitude is shown in 
units of the satellite radius (Tab QJ) . Vertical error bars reflect statistical uncertainty due to the small 
number of impacts. The straight lines are mean measured curves for each moon (Table IIVI col. 8) 
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Table I: Physical properties of Jupiter and the Galilean satellites used in this paper. The Hill radius 
is defined as rnm = '"(ap^^)^''^ with M, m being the masses of Jupiter and the moon, separated by 
distance r. 



Object 


Jovicentric 


Radius 


Symbol 


Hill radius 


Escape speed 




distance 


^obj 




^Hill 


'^esc 




(Rj) 


(km) 




(Robj) 


(kms^-*-) 


Jupiter 




71,492 


Rj 






lo 


5.9 


1,818 


Ri 


5.8 


2.56 


Europa 


9.4 


1,560 


Re 


8.7 


2.03 


Ganymede 


15.0 


2,634 


Rg 


12.0 


2.74 


Callisto 


26.3 


2,409 


Rc 


20.9 


2.44 
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Table II: Galileo flyby characteristics and parameters for the dust particles detected during Galileo's 
flybys at Europa, Callisto and lo (for Europa, particles within 8 Re altitude have been included, and 
for Callisto and lo within 6Rc and 6Ri, respectively): Flyby number (col. 1), time of flyby (col. 2), 
altitude at closest approach to satellite (col. 3), velocity of Galileo relative to s atellite (col. 4), average 
measured particle velocity (velocity error factor VEF < 6, iGriin et al.1 ()l99,^ : col. 5), spin-averaged 
sensor area (maximum value 235 cm^; col. 6), number of class 2 and class 3 satellite particles for which 
their complete data set has been transmitted to Earth (180° < G < 360°; col. 7), number of all events 
(dust plus noise) detected with 180° < © < 360° (both within 8 Re, col. 8), completeness of data 
set due to incomplete data transmission of Galileo (col. 9) and slope of power law fit to the radial 
variation of the impact rate (Fig. 0J 121 1131 col. 10). The slopes are weighted with the square root 
of the number of particles, and the uncertainty takes into account the error bar of each data point. 
Values given in parentheses have been derived from only four or fewer particles. 
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(km) 
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(%) 
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(7) 


(8) 


(9) 


(10) 












Europa 










E4 


96-354.287 


698 


5.7 


6.8 ±4.0 


233 


18 


23 


86 


-2.21 ±0.36 


E6 


97-051.713 


586 


5.7 


5.9 ±4.0 


224 


3 


7 


88 


(-1.69 ±0.58) 


Ell 


97-310.855 


2,043 


5.6 


5.5 ±2.2 


235 


10 


11 


90 


-1.43 ±2.25 


E12 


97-350.502 


201 


6.3 


4.5 ±1.7 


142 


12 


13 


61 


-2.72 ±0.91 


E14 


98-088.556 


1,644 


6.5 


3.4 ± 1.3 


61 


3 


4 


100 


(-1.26 ±0.44) 


E15 


98-151.884 


2,515 


6.4 


5.5 ±1.5 


155 


4 


6 


100 


(-2.11 ±0.59) 


E17 


98-269.163 


3,582 


6.0 


2.8 ±0.8 


180 


4 


4 


54 


(±0.27 ± 1.56) 


E19 


99-032.097 


1,439 


5.8 


5.7 ±4.7 


152 


lot 


13+ 


100+ 


-1.70 ±0.62 












Callisto 










C3 


96-309.566 


1,118 


8.0 


6.3 ±1.7 


224 


10 


35 


80 


-1.90 ±0.45 


C9 


97-176.575 


415 


8.2 


7.6 ±0.2 


197 


3 


5 


100 


(-0.69 ±0.10) 


CIO 


97-260.013 


538 


8.2 


6.3 ±2.3 


142 


22 


37 


64 


-1.31 ± 1.06 


10 


95-341.740 


892 


15.0 


10.3 ±8.4 


lo 

133 


4 


13 


87 


(-2.00 ±1.88) 



f: Data transmission ceased 1 Feb 99, 02:38:46h, i.e. 19 min after closest approach, at 6,767km 
altitude. All data sets of particles detected earlier during the E19 encounter were transmitted. 



Table III: Model parameters for different satellites: speed of impactors, geometric albedo, assumed 
silicate content, energy partitioning parameter, characteristic yield, parameters of the ejecta speed 
distribution (see Krivov et al. 2002b for description of the parameters). 



Satellite 




A 


Gsii 




Y 


no 


7 




( kms""*^) 






(%) 




( m s^"*") 




lo 


26 


0.61 





30 


2.8 X 10^ 


28 


2.0 


Europa 


21 


0.64 





30 


1.6 X 10^ 


30 


2.0 


Callisto 


15 


0.20 


70 


20 


7.1 X 10^ 


51 


1.4 
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Table IV: D e rived properties of the impact-generated dust clouds (data for Ganymede taken from 
Kriiger et al. (|2f)flfl h): mean flyby speed of Galileo at the satellite (col. 2), number of detected cloud 



particles for which the complete data set has been transmitted (col. 3), mean measured impact speed of 
these particles (VEF < 6; col. 4), mean particle mass (col. 5), corresponding particle radius assuming 
spherical particles with density Igcm^'^ (col. 6), slope of the power law fit to the cumulative mass 
distribution (col. 7) and average slope of radial number density distribution (col. 8; averages of slopes 
from individual flybys). The speed of Galileo relative to the satellite has been assumed for the values 
in col. 5 to 7. For the uncertainties in the number density slopes (col. 8) the uncertainties of the slopes 
from each individual flyby have been taken into account. 



Object 


Average 


Number 


Average 


Average 


Average 


Slope a 


Slope of 




flyby 


of 


particle 


particle 


particle 


of mass 


radial number 




speed 


detections 


speed 


mass 


radius 


distrib. 


density distrib. 




(kms~^) 




(kms^^) 


(kg) 


(^m) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


lo 


15.0 


4 


10.3 ±8.4 


8.5 X 10^1" 


« 0.6 




(-2.00) 


Europa 


6.0 


64 


5.5 ±3.5 


1.3 X lO"!-* 


« 1.0 


0.58 ±0.04 


-2.02 ±0.63 


Ganymede 


8.2 


38 


7.2 ±4.9 


9.5 X IQ-^^ 


« 0.6 


0.82 ±0.04 


-2.82 ±2.60 


Callisto 


7.2 


35 


6.4±2.1 


5.2 X IQ-^^ 


« 0.5 


0.86 ±0.14 


-1.60 ±0.39 



Table V: Mass budgets of the dust clouds (model estimates, see Krivov et al. 2002b for description 
of the model). 



Satellite 


lo 


Europa 


Ganymede 


Callisto 


Mass flux of impactors [gm~'^s~"^] 


1 X 10"^^ 


7 X 10"^^ 


4 X 10"^^ 


3 X 10"^^ 


Mass inflow of impactors [gs^"*^] 


130 


50 


100 


50 


Yield 


3 X lO'^ 


2 X 10^ 


4 X 10^ 


7 X 102 


Mass production rate of ejecta [gs^-*^] 


4 X 10^ 


9 X 10^ 


4 X 10^ 


3 X 10^ 


Mean lifetime of ejecta [s]") 


50 


70 


120 


250 


Steady-state mass of the cloud [tons] 


200 


60 


40 


9 


Fraction of escaping ejecta 


1 X 10"^ 


2 X 10-'' 


8 X 10""^ 


4 X 10-3 


Ejection rate into circumjovian space [gs^^] 


400 


200 


300 


100 



Residence time within the Hill sphere before recoUision or escape 



46 



